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Method and arrangemoit for optical coherence tomography 



FIELD OF THE INVENTION 

The present invention relates to the field of medical imaging, and more 
particularly to optical coherence tomogr^hy (OCl^. 

BACKGROUND AND PRIOR ART 

OCT is a known technique for non iavasive cross-sectional imaging in 
biological systems. OCT uses coherence interferometry to produce a two-dimensional image 
of optical scattering &om internal tissue micro structures in a way that is analogue to ultra 
sonic pulse-echo imaging. 

In coherence reflectrometry, the coherence property of light reflected from a 
sample provides information on the time-of-flight delay from the reflective boundaries and 
back scattering sites in the sample. The delay information is then used to determine the 
longitudinal location of the reflection sites. The OCT system performs multiple longitudinal 
scans at a series of lateral locations to provide a two-dimensional map of reflection sites in 

the sample. Corresponding OCT systems are known from US patent no. 5,847,827 and US 

patent no. 6,057,920 and J. Biomed. Opt 4(1), 157-173 (1996) . 

A common disadvantage of prior art OCT systems is that a high lateral 

resolution is inevitably associated with a strong decrease in available scanning depth. 

The present invention therefore aims to provide an improved method for optical coherence 

tomography and a corresponding computer program product and arrangement for optical 

coherence tomography. 

SUMMARY OF THE INVENTION 

The present invention provides for a method for optical coherence tomography 
whereby the combined Ught beam which results from the combination of the reference light 
beam and the measurement light beam is sampled. For each of the sampling positions an 
amplitude of the intensity variation due to interference is determined when the reference light 
beam is modulated. The amplitude is also referred to as modulation depth. 
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The ampUtudes of the individual sampling positions are then added which 
provides an intensity signal for one picture element. This intensity signal has an improved 
quality, which is due to the feet that adding ampUtudes while neglecting phase information 
allows the use of a higher collection angle (or numerical aperture), without the reduction in 
depth range that occurs in prior art oct systems with a higher NA. . In particular the effect of 
speckle can be reduced this way (cf. Journal of Biomedical Optics 4(1), 95-105 (January 
1999), "Speckle in Optical Coherence Tomography", J. M. Schmitt, et al.). 

In accordance with a preferred embodiment of the invention the phase 
information of the intensity variations at the sampling positions is used for compensating the 
respective phase offsets. This is done by controlUng an adjustable optical filter having an 
adjustable optical element for each of the sampling positions. 

This way an optical system having a high numerical aperture can be used 
without limiting the scanning depth to the focal region of the optical system, while fully 
exploiting the resolution gain due to the high numerical aperture. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the foUowing preferred embodiments of the invention wiU be described in 

greater detail by making reference to the drawings in which: 

Figure 1 is a block diagram of an embodiment of the present invention. 
Figure 2 is illustrative of an OCT system using a bundle of optical fibres. 
Figure 3 is illustrative of an OCT system using a single optical fibre for 

receiving the combined Ught beam, while sequentially sampling the combined beam. 

Figure 4 is a schematic diagram illustrating the image acquisition for pixels 

outside the focal plane. 

Figure 5 is a block diagram of a signal processing system, 
Figure 6 shows an image acquired in accordance with the principles of the 
present invention in comparison to prior art image acquisition. 

DETAILED DESCRIPTION 

Figure 1 shows an OCT system 100. The OCT system 100 has a hght source 
102. The low-coherence light beam 104 provided by light source 102 is directed on beam 
splitter 106 to provide a reference light beam 108 and measurement light beam 1 10. In 
operation the light source 102 outputs low-coherence light. 
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Reference Ught beam 108 can be modulated by longitudinal movement of 
reference mirror 1 12 and / or by modulator 1 14. The modulation of the reference Ught beam 
108 by means of reference mirror 1 12 and / or by means of modulator 114 is as such known 
from the prior art. 

Alternatively the modulation can be performed by modulating the frequency 
of the reference light beam 108 by means of modulator 1 14. In this instance an acousto-optic 
modulator can be used for providing frequency modulation in the range of 1 kHz or up to 1 
MHz. 

Measurement light beam 1 10 is focussed on point 1 16 in object 1 18 by means 
of scan unit 120 and objective lens 122. The reflected measurement Ught beam 1 10 and the 
reflected reference light beam 108 are superimposed to provide combined Ught beam 124. 
The modulation of the reference beam results in an intensity modulation in the combined 
beam due to interference only when there is a reflecting site inside the imaged object 1 1 8 
such that the difference in optical path length between on the one hand the light that has 
travelled to that site and back to the beam splitter 1 06. and on the other hand the Ught that has 
travelled to the reference mirror 1 12 and back to the beam spUtter 106 is less than the 
coherence length of the Ught source 102. Thus, the reference mirror selects a window mside 
the object with a width equal to the coherence length of the Ught source. 

Combined Ught beam 124 is sampled in a plane which is perpendicular to 
combined light beam 124. An array of photo detector elements 126 is provided for sampUng 
combined Ught beam 124. Each of the photo detector elements of the array 126 deUvers an 
electrical measurement signal for the intensity of the Ught impinging on the surface of the 
corresponding detector element. Each of the electrical measurement signals is transmitted by 
a separate Une 128 to signal processing unit 130. Signal processing unit can be reaUzed as an 
analogue or a digital device. 

Signal processing unit 130 provides an amplitude signal for each one of the 
photo detector elements via corresponding lines 132 to image processing unit 134. 

The corresponding wave front of the measurement light beam 1 10 is focussed 
on point 1 16 and reflected back from object 118. When this reflection occurs the wave front 
profile is altered. This effect is also referred to as aberration. The magnitude and shape of the 
aberration depends on the various light paths through object 1 1 8 of the measurement hght 
beam 1 10. As a consequence the interference of the reflected measurement Ught beam 1 10 
with the reflected reference light beam 108 is non-unifoim due to the modified wave front 
profile of the reflected measurement Ught beam. 
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^^""«^^^^^b«<«°^<^tive interference at a sampling position Of on^ 
photo detector element of array 126 while there is destructive interference at another 
sampling position at the same time. When the reference light beam 1 08 is modulated this 
translates into a corresponding non-unifomx intensity variation at the various sampling 

► positions provided by array 126. 

For each one of the sampling positions an amplitude of the intensity variation 

dunng the modulation ofthe reference Ught beam 108 is measured. All ofthe measured 
amphtudes are added irrespective of the phases of the intensity variations to provide a single 
intensity signal. This intensity signal is the resulting image information for point 1 1 6. 

If the signals fiom Ihe detector elements were added directly (as would be the 

case for one single large area detector as used inprior art systems)intensity variations^ 
position would cancel out intensity variations with opposite phase at anotiier position 
resulting in signal loss and information loss. 

By means of scan unit 120 object 118 can be scamied in order to provide a 
tomographic unage by image processing unit 134. 

fix addition the signal processing unit 130 can also provide phase information 
for tiie mtensity variations of the individual photo detector elements. For tins purpose an 
ari^itrary reference phase can be defined. Ttc. phase offset of each of ttxe intensity variations 
of the photo detector elements of array 126 due to tire modulation of tire reference tight beam 
108 rs provided from signal processing unit 130 to control unit 136. Contiol unit 136 is 
coupled to adaptive optical filter 138. 

Adaptive optical filter 138 has an array of adjustable optical elements. Each 
one of tire adjustirble optical elements is assi^ed to one of tire sampling positions provided 
by array 126. For example tire adjustable optical elements can be realized by means of a 
controliable LCD element. 

Control mut 136 controls the adaptive optical filter 138 in order to compensate 
the phase offiets which are determined by the array 126. 

For example, the reflected measurement Hght beam 1 10 comprises Hght path 
140 going through optical element 142 of adaptive optical filter 138. Light path 140 extends 
to beam spUncr 106 ftom where it goes to photo detector elem«,t 144 of flre array 126 In 
other words hght path 140 defines a one to one relationship between opdcal element .42 «,d 
photo detector element 144. The Ught which impinges on photo detector element 144 via 
hght path 140 is measured by photo detector element 144 in order to detennine the local 
amphhide of tire interference pattern. 
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The phase offset of the intensity variation of the measurement signal provided 
by photo detector element 144 is provided to control unit 136. Control unit 136 has a suitable 
control filter such as a proportional integral (PI) or a proportional integral differential (PID) 
filter to provide a compensation phase. A corresponding control signal is outputted by control 
unit 136 such that the optical element 142 is contirolled to apply this compensation phase to 
the Ught of measurement light beam 110 going through optical element 142 on light path 140. 
This way a control loop is formed in order to compensate the change of the wave front profile 
which is caused by the reflection of the measurement Ught beam 1 10 on the object 1 1 8. 

It is a particular advantage of using a control loop for compensating the 
modification of the wave front profile that aberration is thereby corrected and defocus is also 
compensated. This enables to use an objective lens 122 having a large numerical aperture for 
improved resolution. The geometiical focus will in this case automatically follow the depth 
scan provided by the reference mirror 1 12. 

It is to be noted that defonnations and aberrations can be corrected either in 
the electronic, or in tiie optical donuiin. 

Figure 2 shows an alternative embodiment. A fibre bundle 200 is located 
inside the image of the objective lens pupil 202. The fibre bundle 200 can be used as an 
adaptive optical filter by selectively applying stress on individual fibres of the fibre bundle 
200. By application of stress onto an individual fibre of the fibre bundle 200 a phase shift of 
the Ught going through that fibre is induced. Hence, fibre bundle 200 has a fimction which is 
analogous to optical element 142 of figure 1. For simplicity the Ught source, the reference 
light beam and the combined Ugjit beam with the associated optica] elements are not shown 
in the embodiment of figure 2. 

Figure 3 is illustrative of a sequential sampUng set-up. 
As opposed to the embodiments of figures 1 and 2 the sampling is not 
performed in parallel by using multiple photo detector elements but sequentially using a 
single optical fibre 300 as a receiver for the combined light beam (c£ combined Ught beam 
124 of figure 1). In the embodiment considered here the sampling is performed by pivoting of 
scamring mirror 302. This way the scanning muror moves the limited numerical aperture of 
the OCT set-up across the higher numerical aperture of the objective lens 304. In the case of 
a moving object (in-vivo measurements) motion compensation algorithms can be used to 
correct for artefacts that occur when the object is moved during the scan over the combined 
beam. 
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Figure 4 is illustrative of the principles for generating an image for image 
plane 400 which is displaced from focal plane 402 of the objective lens by a depth d. 

In order to obtain an intensity value for the picture element x' on image plane 
400 a lateral scan is performed on focal plane 402. In particular the lateral scan is performed 
in order to obtain intensity signals for the picture elements x,, xj, . . . Each of the picture 
elements x; has an assigned Hght cone containing the picture element x' on image plane 400. 
For example picture element x, on focal plane 402 has Ught cone 404 containing the picture 
element x'; likewise picture element xa on focal plane 402 has Ught cone 406 which also 
contains the picture element x'. 

Light cone 404 contains reflected light which is detected on a position y, when 
the sampling of the combined Ught beam for acquiring an intensity signal for the picture 
element xi is performed. Likewise designates another sampUng position for the 
measurement of Ught cone 406 when the data acquisition to provide the intensity signal for 
the picture element Xz is performed. In general each one of tiie pichire elements Xj on the 
focal plane 402 has a corresponding sampling position yj for measurement of the Ught cone 
containing the picture element x'. 

If an array of detector elements is used as it is the case in the embodiments 
figures 1 and 2 y, designates one of the photo detector elements of the array (cf. array 126 
and photo detector element 144 of figure 1). 

The intensity signal for tiie picture element x' is obtained by adding of the 
ampUtudes obtained for the light cones 404, 406 .. . during the sequential lateral scan within 
focal plane 402. By doing tiiis for all picture elements on image plane 400 this way an image 
can be calculated for image plane 400 which is distanced from the focal plane 402 even 
without using the adaptive optical filter of figures 1 and 2 (adaptive optical filter 138 of 
figure 1 and fibre bundle 200 of figure 2). 

Expressed in mathematical terms the intensity signal I for the picture element 
x' is obtained as follows: 

I(x',d) = Sumj [ Sj ( x'-(d/f)yj . d) ] 

where 



I(x',d): Intensity signal for picture element x' 

j • Index for sampUng position 
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yj ■ Location of the jfli sampling positon 

Sj(xi,d): Intensity signal if only data fiom the jth position is used 

Foc^ length of the objective lens 
X*: x'=xi + (d/f)xi 

Figure 5 is illustrative of an electronic device 500 for optical coherence 
tomography in accordance with the present invention. Electronic device 500 can be realized 
by a number of separate physical units or it can also be realized as a "system on chip". 

Electironic device 500 has an array of photo detector elements 502 (cf. array 
126 of figure 1). Each of the array elements yj receives a portion of the combined Ught beam 
on a corresponding light path j. 

When the reference Ught beam is modulated the ampUtude of the intensity 
variation is measured for each of the sampling positions y, . The amplitude information of all 
the sampling positions is added by integrator 504. The resulting intensity data is outputted as 
an image data for the actual picture element and stored in storage 506. 

Further the phase offset of each of the intensity variations on the sampling 
positions y, is inputted into phase control 508. If an adaptive optical filter is available in the 
light path of tiie measurement light beam, phase. controlAOS can.form a control loop with the 
adaptive optical filter in accordance with the principles as explained with respect to figure 1 
in order to compensate aberration and defocus. 

If such an adaptive optical filter (cf. adaptive optical filter 138 of figures 1) is 
not available the approach as explained with respect to figure 4 is used as an alternative for 
calculating image data for image planes which are outside the focal plane. For this purpose 
all of the amplitudes of the intensity variations measured at the sample position y, are stored 
with the corresponding pichne element X| as a key in storage 510. By means of this data basis 
integrator 504 evaluates the above formula 1 in order to calculate intensity data for picture 
elements x' which are on an image plane outside the focal plane. Again this image data is 
stored in storage 506. 

Figure 6 illustrates the beneficial effect of flie present invention. 

The image a) of figure 6 has been acquired by means of OCT without 
compounding. Image b) has been obtained by averaging different images but also without 
compounding. Image c) has been obtained by averaging across different regions with 
compounding in accordance with the embodiment of figure 3 without usage of an adaptive 
optical filter. Speckle noise in the image c) is clearly suppressed in comparison with unages 
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a) and b). In particular, the image b) demonstrates that averaging of images alone is not 
sufficient to suppress speckle. 
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100 OCT system 

102 light source 

1 04 low-coherence light beam 

106 beam spUtter 

10 108 reference Ught beam 

1 1 0 measurement light beam 

112 reference mirror 

114 modulator 
116 point 
15 118 object 
120 scan imit 

1 22 obj ective lens 

1 24 combined light beam 

126 array of photo detector elements 

20 128 line 

130 signal processing unit 

132 line 

1 34 image processing unit 

136 control unit 

25 138 adaptive optical filter 

140 light path 

142 optical element 

144 photo detector element 

200 fibre bundle 

30 202 objective lens pupil 

300 optical fibre 

302 scanning mirror 

304 objective lens 

400 image plane 
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402 


focal plane 


404 


Ught cone 


406 


light cone 


500 


electronic device 


502 


array 


504 


integrator 


506 


storage 


508 


phase control 


510 


storage 



